Anomalous isotope effect in phonon-dressed iron-based superconductors 
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The role of electron-phonon interactions in iron-based superconductor is currently under debate 
with conflicting experimental reports on the isotope effect. To address this important issue, we 
employ the renormalization-group method to investigate the competition between electron-electron 
and electron-phonon interactions in these materials. The renormalization-group analysis shows 
that the ground state is a phonon-dressed unconventional superconductor: the dominant electronic 
interactions account for pairing mechanism while electron-phonon interactions are subdominant. 
Because of the phonon dressing, the isotope effect of the critical temperature can be normal or 
reversed, depending on whether the retarded intra- or inter-band interactions are altered upon 
isotope substitutions. The connection between the anomalous isotope effect and the unconventional 
pairing symmetry is discussed at the end. 



What role phonons play for electron pairing always 
stimulates intense discussions when a new supercon- 
ductor is discovered [1 such as cuprates[2 , magnesium 
diboride[3 , sodium cobalt oxide [4] and recently found 
iron-based superconductors [51 [6]. The first checking 
point is the critical temperature of superconductivity. 
By studying density of states for phonons, the critical 
temperature can be estimated from the Bardeen-Cooper- 
Schrieffer (BCS) theory or the Migdal-Eliashberg formula 
and then compared with the experimental data to see the 
deviation from conventional superconductivity[7-ll . In 
unconventional superconductors, the phonon-mediated 
interactions are insufficient to explain the pairing mech- 
anism and it is of crucial importance to study the in- 
terplay between electron-electron and electron-phonon 
interactions [TJHlT]. For instance, even when the pairing 
mechanism is electronic origin, dispersions observed in 
angle-resolved photoemission spectroscopy manifest dis- 
tortions upon isotope substitutions [TBHSl]- 

To reveal the role of phonons, isotope effect in the su- 
perconducting temperature Tc is commonly used[2[[3[|25]- 
[30j. By detecting the change in Tc upon isotope substi- 
tutions, the isotope exponent a = — (i(logTc)/(i(log M) 
can be obtained [STJ [32], where M is the mass of the 
substituted element. In conventional superconductors 
where electron-phonon interactions reign, the exponent 
is a = 1/2. On the other hand, in an unconventional 
superconductor without any relevant phonon-mediated 
interactions, a = is expected. The isotope effect ob- 
served in iron-based superconductor [161 IIZl 1211 I27H3Q] 
seems to tell a more complicated story. For instance, a 
strong isotope effect by iron substitution [27] is found in 
SmFeAsOi_a;Fa; and Bai_a;KajFe2As2, almost as large as 
that in conventional superconductors. On the contrary, 
inverse isotope effect [28] is also spotted in (Ba,K)Fe2As2 
with different isotope substitutions. 

Motivated by the controversy, we investigate the com- 



petition between electron-electron and electron-phonon 
interactions by the unbiased renormalization-group (RG) 
method. Due to the retarded nature of the phonon- 
mediated interactions, the energy dependence must be 
included. The minimal approach to include both simul- 
taneous and retarded interactions can be accomplished 
by the step-shape approximation [33l l34] as shown in Fig. 
1(a), 

9i{^) = gSi^D - (1) 

where gi and gi represent (instantaneous) electronic in- 
teractions and (retarded) phonon-mediated ones. The 
energy scale for the retarded interactions is set by the 
Debye frequency ujd- Our RG analysis reveals that the 
pairing mechanism is dominated by the electronic inter- 
actions gi. But, the retarded interactions gi also grow 
under RG transformation and become relevant in low- 
energy limit. Inclusion of these subdominant interactions 
leads to anomalous isotope effect. Though a quantitative 
estimate of the isotope exponent a is difficult in the pres- 
ence of both types of interactions, we demonstrate how 
the exponent can be extracted numerically from RG flows 
in weak coupling. Surprisingly, the sign of the exponent 
a sensitively depends on whether the inter- and/or intra- 
band interactions are altered by isotope substitutions. 

To illustrate how the RG scheme works, we start with 
a five-orbital tight-binding model for iron-based super- 
conductors with generalized on-site interactions, 

p,ci,b OL i^a 

+ ^2 X ^riiaariib/S ^ Jh X y^.4aa^iba4^0Ciai3 

i,a<b a, (3 i,a<b a, (3 

^ JhYI [4at4aiCinCi6t + H.c] I, (2) 
i,a<b ^ 

where a^b = 1,2, ...,5 label the five d-orbitals of Fe, 
1 : dsz'^-R^, 2 : dxz, 3 : dyz, 4 : dx^-v'^, 5 : dxY, 
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FIG. 1: (a) Step-like interaction profile for simultaneous and 
retarded interactions. A sharp step is assumed at the Debye 
frequency ujd- (b) Fermiology of the five-band model x = 0.1. 
These Fermi surfaces are well sampled by five pairs of Fermi 
points, equivalent to the four- leg geometry with quantized 
momenta (dashed lines). 
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FIG. 2: RG exponents for (a) the simultaneous and (b) the 
retarded Cooper scatterings. The dominant interactions are 
pairing hopping between and within band 1 and band 2, with 
maximal exponent of unity, while other relevant couplings are 
subdominant with RG exponent smaller than one. 



and a =t, i is the spin index. The kinetic matrix Kab 
in the momentum space has been constructed in previ- 
ous studies [3F. The generalized on-site interactions con- 
sist of three parts: intra-orbital /7i, inter-orbital U2 and 
Hund's coupling Jh- Adopted from previous studies, we 
choose the values, Ui = 4 eV, U2 = 2 eV and Jh = 0.7 
eV for numerical studies here. 

Fermiology is important in the multi-band supercon- 
ductors. The electron doping x is related to the band 
filling n = 6 + X (n = 10 for completely filled bands) 
here and the Fermi surface aix = 0.1 is illustrated in 
Fig. 1(b). There are five active bands: two hole pock- 
ets centered at (0,0) and another hole pocket centered 
at (7r,7r) while two electron pockets located at (7r,0)and 
(0,7r) points [36 . To simplify the RG analysis, we sample 
each pocket with one pair of Fermi points (required by 
time-reversal symmetry). This is equivalent to a four- 
leg ladder geometry with quantized momenta as shown 
in Fig. 1(b). In the low-energy limit, the effective 
Hamiltonian j37] - [39] is captured by five pairs of chiral 
fermions with different velocities. 

The interactions between these chiral fermions fall into 



two categories [341 [4T] : Cooper scattering c-^ , and for- 
ward scattering flj , ffj . The retarded ones share the 
same classification, denoted with an extra tilde symbol. 
The RG equations for the simultaneous interactions are, 
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k^kj ^ ^ik^kjj 
2 



(4) , 



2r- r^. 



(3) 



where g = dg/dl^ where / = ln(Ao/A) is the logarithm of 
the ratio between bare energy cutoff Aq and the running 
cutoff A. The tensor aij^k = {^i + Vk){vj + Vk)/[2vk{vi + 
Vj)] with Vi representing the Fermi velocities. 

The second set of equations describes how the retarded 
interactions are renormalized. 
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Note that we separate the intra-band and inter-band cou- 
plings for clarity, i.e. i ^ j in the above RG equations. 
In fact, the separation is necessary because we shall see 
later that inter-band and intra-band couplings play dif- 
ferent roles in the low-energy limit. In addition, fa = 
and fa = to avoid double counting. 

By integrating the two sets of RG equations numeri- 
cally, we found all couplings are well described the scaling 
ansatz |42] , 



Gi 



{Id - iV^^ ' 



(Id - iv~^^ 



(5) 



where Gi^Gi are non-universal constants and ^g^-,^g^ are 
RG exponents for simultaneous and retarded couplings. 
The divergent length scale Id^ associated with the pairing 
gap, is solely determined by electronic origin. The dom- 
inant pairing occur within band 1 and band 2 and the 
Cooper scatterings cn, C22 C12 have maximum exponent 
7i = 1. Other Cooper scatterings are subdominant with 
exponents close to 0.9, as shown in Fig. 2(a). Mean- 
while, by Abelian bosonization[38| [39J, the signs of 
from numerics lead to sign-revised (between electron and 
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FIG. 3: (a) Interaction profile of the dominant intra-band 
C and inter-band C± Cooper scatterings before RG trans- 
formation (b) As RG progresses, the step evolves since the 
Debye energy is rescaled, ujd{1) — ujdg^ - (c) For / > Id, 
the distinction between simultaneous and retarded interac- 
tions disappears, (d) Schematic picture for phonon-dressed 
unconventional superconductor. 



hole pockets) 5-j--wave pairing, agreeing with the previ- 
ous functional RG study [40]. Note that these exponents 
are rather robust within the doping range where the same 
Fermiology maintains. What about the phonon- mediated 
interactions? As clearly indicated in Fig. 2(b), the 
RG exponents for Cii,C22 are roughly 0.6, much smaller 
than the dominant electronic interactions, showing the 
pairing mechanism is electronic origin. However, since 
the RG exponents are positive, the retarded interactions 
also grow under RG transformation. These subdominant 
phonon-mediated interactions can lead to anomalous iso- 
tope effect as explained in the following. 

To achieve quantitative understanding in weak cou- 
pling, the rescaled Debye frequency must be taken into 
account carefully. Under RG transformations, ud 
uod^^ as shown in Fig. 3. At the (logarithmic) length 
scale Id = \og{Ao /ujjj)^ the difference between gi and 
gi disappears. The Debye frequency cji) ~ 30 meV in 
iron-based materials [11 and the band width (thus Aq) 
is 3 — 4 eV, giving rise to Id ^ 5. Note that the RG is 
truncated at the cutoff length scale Ic where the maximal 
coupling reaches order one. In weak coupling, it is clear 
that Ic > Id and thus the RG scheme must be divided 
into two steps. For I <Id^ both sets of RG equations are 
employed. At I = Id ^ the functional form for the retarded 
interactions is the same as the instantaneous one. Thus, 
one should add up both types of couplings gi{lD)-\-gi{^D) 
and keep running RG by just the first set of equations. 
In physics terms, this means that the difference between 




FIG. 4: The cutoff length scale 
intraband interactions ^^(O) = - 



versus Id for inclusion of 
-0.3[/ (blue circles) and in- 



terband ones Cij{^) — — 0.14[/ (red square), where U is the 
strength of electron-electron interactions. For convenience, 
the axes are rescaled in the unit of /e, the cutoff length scale 
with electronic interactions only. The inset shows the isotope 
exponent by taking numerical derivative. 



simultaneous and retarded interactions vanishes before 
the pairing gaps open. 

Numerical results for the two-step RG indicate the 
same superconducting phase as described in previous 
paragraphs but the isotope exponent a can be ex- 
tracted numerically. Since ujd ^ M~^/^, it means that 
(i(logM) = —2d{\ogujD) = 2dlD' By standard scal- 
ing argument, the critical temperature takes the form, 
ksTc ~ A [^(0)] = Ace~^^, where Ac is the pairing gap 
at the cutoff length scale. With straightforward algebra, 
the isotope exponent can be written as 



a = 



^(logTc) _ I d(logAc) 1 die 
~d{\ogM) ^ 2d{\ogujD) 2dlD 



(6) 



For conventional superconductor, /Ac ^ and the cut- 
off length scale is not sensitive to the Debye frequency 
(the second term vanishes). Thus, a ^ 1/2. On the 
other hand, for unconventional superconductors without 
relevant electron-phonon interactions, Ac ~ Aq and the 
cutoff length scale is also not sensitive to the Debye fre- 
quency. It is clear that a = in this case. 

But, what happens if the electron-phonon interactions, 
though not dominant, are actually relevant under RG 
transformation? We study how the cutoff length scale 
Ic varies with different Debye frequencies due to isotope 
substitutions. In weak coupling, we found that g^ are 
much larger than gi. Thus, Ac has very weak dependence 
on uod and the first term can be ignored. The contribu- 
tion from the second term is shown in Fig. 4. We tried 
two different profiles for the retarded interactions. In- 
clude only intra-band interactions, Qi(0) = —0.3/7 first, 
where U is the strength of electron-electron interactions. 
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The isotope exponent is positive (reading from the slope) , 
a ~ 0.1, with very smooth variation. On the other hand, 
with only inter-band interactions, Qi(0) = —0.14/7, the 
isotope exponent is negative and changes gradually from 
zero to a ~ —0.03. These anomalous isotope effect is 
closely related to the pairing symmetry. For the 5±-wave 
pairing, ca < but Cij > at the cutoff length scale. The 
phonon-mediated intra-band interactions ca < help to 
develop the pairing instability and thus lead to a posi- 
tive isotope exponent. On the other hand, the inter-band 
ones Cij < have opposite sign with their simultaneous 
counterparts Cij. In consequence, the pairing instability 
is suppressed and an inverses isotope effect is in order. 
The RG analysis presented here provide clear and natu- 
ral connection between the anomalous isotope effect and 
the unconventional pairing symmetry. 

Although the isotope exponent a can be extracted 
numerically in weak coupling, extending the quantita- 
tive description to intermediate coupling may not be 
easy. If the pairing gaps open before hitting the De- 
bye energy scale, i.e. Ic < Id^ our numerical results 
show that Ic solely depend on electronic interactions 
and thus dlc/dln = 0. The isotope exponent in this 
regime mainly arises from the first term. The pairing gap 
Ac = Ac(Ao, ct;De^), depending on both the bandwidth 
and the rescaled Debye frequency, is now quite compli- 
cated. The RG analysis alone is not sufficient to obtain 
a in a quantitative fashion. However, we recently found 
that the effective Hamiltonian at the cutoff length scale is 
well captured by mean- field theory [43 . In principle, one 
can combine RG and mean-field approaches together to 
compute the isotope exponent in intermediate coupling. 

In summary, we investigate the competition between 
electronic and phonon-mediated interactions in iron-base 
superconductor by RG analysis. The pairing mecha- 
nism (thus pairing symmetry as well) is determined solely 
by electronic interactions. However, we show that the 
phonon-mediated retarded interactions are also relevant 
and give rise to anomalous isotope effect. Due to the 
unconventional s± pairing symmetry, the isotope expo- 
nent sensitively depends on the profile of the retarded 
interactions. 
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